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ABSTRACT 

One of the design requirements of the International Space Station (ISS) is that each year 
accelerations of one micro-g cannot be exceeded at the ISS internal payload location for 6 periods 
of not less than 30 consecutive days. Although there are other causes, this study deals only with 
the accelerations caused by atmospheric drag. The critical ambient neutral density, computed 
using the Marshall Engineering Thermosphere Model, required to produce accelerations of one 
micro-g on the ISS, is estimated using an atmospheric drag acceleration equation. Results show 
that the design requirements may be difficult to meet during periods of extremely high solar 
activity; the planned reboost and altitude strategies for the ISS may have to be revised to allow 
for the uncertainty in the prediction of neutral atmospheric density within the 100-day period 
established for orbital decay before reboost 
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NOMENCLATURE 

F10.7 daily value of 10.7 cm solar radio noise flux, units are 10* 22 watts/m 2 /Hz bandwidth or 
10 4 Jansky 

F 10 7 average value of F10.7 
i inclination, 51.6 deg 

K degrees Kelvin 

r 0 Earth's equatorial radius, 6,378,140 m 
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V c satellite circular velocity, m/s 

Z altitude above the Earth, m 

p. Earth's gravitational parameter, 3.986012 x 10 14 m 3 s' 2 
p neutral density, kg/m 3 

p z density at altitude Z 

co Earth's rotation rate (7.2921 x 10‘ 5 rad/s) 

INTRODUCTION 

Both the micro-gravity experiment acceleration requirements and the altitude and reboost 
strategies of the ISS are dependent on the combination of the ballistic coefficient of the ISS and 
the neutral orbital altitude mass density. Predicting the orbital altitude neutral mass density for 
the time period of the assembly and operation of the International Space Station (ISS) with any 
accuracy is extremely difficult because the exact relationship between the solar inputs that affect 
the atmosphere and the magnitude of the change they cause are poorly known. To make matters 
worse predictions of the proxy solar inputs used in the Marshall Engineering Thermosphere 
(MET) model are even more difficult. 

This study estimates what might occur in the future by simulating ISS operations during 
selected periods in the historical data base of required inputs to the neutral atmosphere model 
used in the ISS development program. Results are compared to the ISS design requirements and 
planned reboost and altitude strategies. The data base of inputs required for the MET model of 
the neutral thermosphere extends from 1947 through 1994. 

BACKGROUND 
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NASA has specified that the MET model of the neutral atmosphere is to be used in the design, 
development, and testing phases of the ISS 1 . The MET model is a computer code 2,3 based on the 
Jacchia 70 4 and the Jacchia 71 5 (for seasonal variation) thermosphere models. The Jacchia 70 4 
atmospheric model documentation states that the density output of the model was referenced to 
the 81 -day mean values of the 10.7 cm solar radio noise flux, F 10.7 • Discussions and studies 
conducted in the middle to late 1970s and a recent study 6 established that use of the 81-day 
running mean prior to the time of application produced predicted lifetimes for Earth orbiting 
satellites that were in better agreement with the actual lifetimes than the use of 162-day average 
10.7-cm values required for implementation of the original version of the MET model 2 . Although 
the decay histories of satellites can be acceptably duplicated after the fact by using the MET 
model with the actual parameters designated in the original documentation, decay histories of 
satellites at future times depend entirely upon the accuracy of (1) the ballistic coefficient and (2) 
the proxy solar parameters required as inputs to the model. Results of studies that were 
conducted to determine what F10.7 values could be predicted for use in the prediction of decay 
rates at future times showed that values of the 10.7 cm solar radio noise flux and the 3 -hourly 
average geomagnetic index averaged over less than 13-months could not be predicted with 
sufficient accuracy for the application. 

DISCUSSION 

The neutral atmosphere parameter that affects the micro-gravity acceleration design requirement 
and the altitude strategy plan of the ISS is the ambient mass density. This analysis of these two 
portions of the ISS program uses the 1947-1994 historical database of the parameters that are 
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required inputs to the MET model. It is in essence a look at the ISS as if it were in orbit at 
selected time periods between 1947 and 1994, primarily periods with high levels of solar activity. 

Calculation of ambient mass density with the MET model is based on determination of the 
temperature structure above 90 km altitude, which is determined by construction of a temperature 
profile that begins at 90 km with a fixed boundary condition (defined below) and ends at about 
350 km where the temperature is defined as the exospheric temperature. The exospheric 
temperature is the key variable for determination of the variability of ambient orbital density. The 
exospheric temperature is highly dependent on the extreme ultraviolet (EUV) output of the Sun 
and the solar storms which inject energetic particles into the Earth’s atmosphere. Since neither of 
these parameters are measurable at the Earth’s surface, proxy parameters that are available are 
used as inputs to the model. The 10.7-cm solar radio noise flux, both the average and the daily 
value, is a proxy for the EUV while the 3 -hourly average geomagnetic index, ap, is a proxy for the 
Joule heating of the atmosphere that results from the deposition of the energetic particles. 

The MET model is used to compute total mass density from 250 to 500 km altitude using the 
available 47 years (1947-1994) daily values for the 10.7-cm solar flux and the 3-hourly (8 per 
day) average values of the geomagnetic index, ap. These mass densities are then used (1) to 
compute micro-gravity levels as a function of altitude for comparison with the specified design 
requirements and (2) in orbital decay analyses for comparison with the altitude strategy plan. 

Once the variations in the exospheric temperature over the globe are computed, the only other 
parameter needed to describe the thermosphere density is the fixed boundary condition at 90 km 
altitude, which is a density value of 3.46 x lO' 6 kg m' 3 , an ambient temperature of 183 K and a 
mean molecular mass M 90 = 28.878 g mole" 1 . 
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The appendix. Table A, gives mass density versus altitude [250 (510) 500 km] and exospheric 
temperature [600 (5200) 2200 K], This convenient reference table for density was computed 
using the MET model for the assigned values of Te. Table A covers the planned operating 
altitude range for the ISS. The values for density match those from the Jacchia 70 model 4 for Te 
greater than 700K; at 600K the differences in density are less than one percent at 500 km altitude. 
This density difference is attributed to the difference in the numerical integration techniques used 
to evaluate the diffusion equation. The density from Table A is illustrated in Fig. 1 . 

(Fig- 1) 

Density variations in the MET model are based on: 

1. The 10.7 cm solar radio noise flux which was measured at Ottawa, Canada from 1947 to June 
1991, and since June 1991 is measured at Penticton BC, Canada. 

2. The planetary geomagnetic index. This index is made available through the International 
Union of Geophysics and Geodesy (IUGG), Gottingen, Germany. 

3. The rotation of the Earth (diurnal variation). 

4. The Earth's rotation about the Sun (semi-annual variation). 

5. Seasonal-latitudinal variations below 170 km altitude. This variation is established by a set 
of empirical equations that correct the mass density between 90 and 170 km altitude. This 
variation does not affect density above 170 km. 

6. Seasonal-latitudinal variations in helium for altitudes greater than 500 km. 

The required input parameters for this application of the MET model are: 

1. Date (month, day, year), 

2. Greenwich mean time (GMT) from midnight, 0000 GMT, in hours and minutes, 
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3. Latitude (deg), 

4. Longitude (deg), 

5. Altitude (km) 

6. The daily value for the 10.7-cm solar flux, Fjo. 7, for the day prior to the date of interest, 

7. The 81-day average 10.7-cm solar flux, fi 0>7 , (81 days prior to day of interest) 

8. The 3-hourly planetary geomagnetic index, ap , 6.7 hours prior to the time of interest. 

The key equations in the MET model that relate to the solar flux and the planetary geomagnetic 
index are presented below. The nighttime minimum global exospheric temperature, Tc, when the 
geomagnetic index, ap, is zero, is: 

T c = 383 + 3.32 Fjo/7 + 1.8(F 107 -F 10 . 7 ) (1) 

where, F10.7 is the daily value and Fi 07 is the average value of the solar flux. In this equation the 
empirical coefficients ensure correct units. 

The daytime maximum exospheric temperature, T D , is: 

T d = RT C ( 2) 

where, R, a function of the 400 day average solar flux, varies from 0.27 to 0.40 and has an 
average value of 0.3 1 4 . 

The variation in the exospheric temperature due to the 3-hourly geomagnetic index, ap, is: 

ATg = 1 .0 ap + 100[1 - exp (-0.08 ap )] (3) 

The ap index has only discrete values between zero and 400 while ATg ranges from 0 to 500K. 
The variable, ATg, is the single most important contributor to exospheric temperature change 
within a day. Three hours is the highest time resolution obtainable from the MET model due to 
the 3-hour intervals for ap. 
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The semi-annual contribution to the exospheric temperature (K) variation is: 


ATs- 2.41 + F 10 7 g(t) 


( 4 ) 


where, 


g(t) = [0.349 + 0.206 sin(360° x + 226.5° )][sin(720° x + 247.6°] 


and 


x = -+.1145| 

y 


l + sin^OP (d/y) + 342 




2.16 


- 1/2 


where, d is days since January 1 and y is the length of the tropical year, 365.2422 days. Fig. 2 
presents the evaluations for g(t). 

(Fig- 2) 

The global maximum exospheric temperature is: 

Tmax = ( 1 . 0+R) T c + ATg + ATs. (5) 

This value occurs at the latitude of the subsolar point with a time lag of approximately 2 hours or, 
at 1400 hours local Sun time. At any altitude in the thermosphere the largest density occurs at the 
location of Tmax. 

The minimum exospheric temperature is: 

Tmin = Tc + ATg + ATs. (6) 

The global averages for density and exospheric temperature in simplified satellite orbital decay 
models are estimated from: 

Po = ~[p(Tmax) + p(Tmin)] (7a) 

Tavg = ( 1 +R/2) T c + ATs + ATg (7b) 
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For applications requiring determination of the density along a satellite trajectory in orbit the 
exospheric temperature as a function of latitude, longitude, and time is required. This local 
exospheric temperature is: 


T L =T C (1+Rsin m e) 


1 + R 


^cos m r| - sin m 0^ 


1 - Rsin m 0 


cos" x / 2 


( 8 ) 


where, x = H + p + psin(H + y) , (-7t < x < 7t) , H is the hour angle of the Sun, m = 2.5, n = 3 .0, 

0 = -37 deg, p = +6 deg, y = +43 deg, and rj = !4 \ LAT- DL| , 0 = 14 | LAT+DL| , 
where, LAT is the latitude of the satellite and DL is the Sun declination angle. 

The local exospheric temperature, Te, is then: 

T E = T L + ATg + ATs (9) 

There are sets of empirical equations that define the temperature versus altitude, T z , between 
the T at 90 km (183K) and Te (at about 350 km) that are functions of altitude and the exospheric 
temperature only. The hydrostatic equation is numerically integrated using these functions, T z , 
and the boundary conditions at 90 km altitude to obtain the neutral mass density up to 105 km 
altitude. For altitudes above 105 km the diffusion equation is used to compute number densities 
of the atmospheric constituents and total mass density. For altitudes above 500 km a correction is 
made to the density for the seasonal variation in helium. 

The value for Te (Eq. 9) lies between Tmin (Eq. 6) and Tmax (Eq. 5), i.e., 

Tmin < Te < Tmax (10). 

The orbit of a low earth orbiting satellite will retrograde approximately 4 degrees per day; 
considering this, there is a good chance that any satellite with an inclination > 23.45 degrees will 
pass through Tmin and Tmax at least one time in 30 days. 
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ANALYSIS 


Fig. 3 presents the daily values for the 10.7 cm solar flux for the period of record under study. 
This period covers the latter part of solar cycle 18, all of cycles 19, 20, 21 and most of cycle 22. 
The 3-hourly planetary geomagnetic index, ap, for this period of record is shown in Fig. 4. The 
daily and 81-day prior average (in this case) solar flux and ap are used to compute the daytime 
daily maximum global exospheric temperature shown in Fig. 5 (top curve); the differences 
between the daily global maximum (Eq. 5) and daily global minimum (Eq. 6) exospheric 
temperature are also shown in Fig. 5 (bottom curve). 

(Figures 3, 4 and 5) 

The largest daily maximum exospheric temperature, 2152.2 K and the largest thermospheric 
density, as computed by the MET model, occurred during solar cycle 19 on September 23, 1957 
at 0600 GMT. The density values versus altitude can be read from Table A for Te = 2150 K. For 
example, at 450 1cm, 400 km and 300 km, the density values are: 0.15699 x 10‘ 10 , 0.26025 x 10' 10 
and 0.83120 x 10' 10 , respectively. Table A can be used to obtain the density for any arbitrarily 
assigned or known value for the exospheric temperature. 

The ISS Microgravity Experiment Requirement from the System Specification Document for 
thelSS 1 is: 

"The Space Station shall provide the following microgravity acceleration performance for 
at least 50 percent of the internal payload locations for 180 days per year in continuous 
time intervals of at least 30 days. At the centers of the internal payloads, a quasi-steady 
(<0.01 Hz) acceleration: (a) magnitude less than or equal to one micro-g and (b) 
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component perpendicular to the orbital average acceleration vector less than or equal to 
0.2 micro-g." 

There are also vibration acceleration limits for the frequency range, 0.01 < f < 300 Hz. 

The variability of atmospheric density during one micro-gravity experiment (30 days duration) 
and the ISS response during this period are important factors that contribute to the residual 
acceleration. 

Atmospheric Drag Acceleration 

The atmospheric drag acceleration, D, (ms' 2 ) is given by: 

D ■ i Cd (£) v *> p (11) 

where Cd is 2.2, A is 2673 m 2 , and m is 420,000 kg for the ISS. These estimates yield: 

C D — = 0.014(m 2 / kg), 
m 

p is atmospheric density and V re i is the relative velocity between the ISS and the atmosphere. For 
a spherical Earth and a circular orbit V re i is a function of the circular velocity (Vc) at altitude (Z), 
inclination (51.6 deg for the ISS) and latitude. At zero deg latitude and an inclination of 5 1 .6 deg, 
V re i at Z = 300 km is 7422 m/s and at Z = 400 km, V re i is 7359 m/s. For illustration purposes V re i 
is taken as 7400 m/s. From the above parameters and assumptions, the atmospheric drag 
acceleration D (Eq. 1 1) can be approximated by: 

D = 383320 p, (12) 

To convert Eq. 12 to micro-gs (p-g) divide by go x 10 " 6 , where go is 9.80665 ms' 2 : 

p-g = 39089 x 10 6 p- (13) 
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For example, when p = 2.5583 x 10' 11 , Eq. 13 gives one micro-g. Hence, when p > 2.5583 x 
10* n , the atmospheric drag acceleration equals or exceeds one micro-g. Fig. 6 shows p = 2.5583 x 
10' n ± 20%, which is a reasonable uncertainty for the MET model, as a function of exospheric 
temperature and altitude. The values in Fig. 6 are important in planning altitude strategies for the 
micro-gravity experiments. The curves in Fig. 6 represent the exospheric temperature to be used 
to estimate the altitude above which one-micro-g will not be exceeded if the MET model 
density : (1) has no uncertainty (curve A); (2) is underestimated by 20 percent (curve B); (3) is 
overestimated by 20 percent (curve C). To protect the ISS for the worst case, which is when 
MET underestimates the density, curve B is applicable. For no uncertainty in the density (curve 
A) one micro-g acceleration would not have been exceeded at altitudes greater than 397 km 
because Te is not greater than 2152.2 K (Fig. 5) for the entire period of record (1947-1994). To 
protect for the worst case uncertainty the altitude would have to be at least 422 km to ensure that 
one micro-g is not exceeded for the entire period of record (curve B). Similarly, if protection is 
required for solar activity levels associated with Te < 1700 K, which excludes very high solar 
activity ,the minimum altitude is 363 km for no uncertainty (curve A) and 380 km for the worst 
case (curve B). 

(Fig. 6) 

A comparison of the mass spectrometer measurements from the Atmospheric Explorer-E (AE- 
E) satellite with the density computed from the MET model using the 162-day centered average 
10.7-cm solar flux, the 81 -day prior average solar flux and the 1/2 gaussian weighted prior 
average solar flux revealed that a better agreement exists between the MET model and the AE-E 
densities using the 81 -day average flux or the 81 -day gaussian average flux than the 162-day 
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average solar flux 6 . The gaussian average flux, the daily flux and ap shown in Fig. 7 for July 1959 
were used in the MET model for the analysis of the ISS in circular orbit at 300 km altitude (Fig. 
8). As shown the density exceeded the critical value for one-micro-g for most of the orbits for the 
entire month of July 1959. This example is for high solar activity. For periods of low solar 
activity one-micro-g is not exceeded for 30 or more days. For this paper T E > 1500K is 
considered high solar activity, T E < 1000K is low solar activity and the values between 1000 and 
1500K are intermediate solar activity. 

(Figures 7 and 8) 

The Variability of Density within An Orbit 

The variability of atmospheric density within an orbit is required for use in ISS engineering 
analyses of (1) the control moment gyros (CMGs’) capability to control the torques produced by 
the difference between the center of pressure and the center of gravity; and (2) the atmospheric 
drag acceleration that may exceed one-micro-g. Fig. 9 for July 15-16, 1959 shows that the 
variation in density within an ISS orbit at 300 km altitude can exceed 2.6 x 10 11 kg m 3 or one- 
micro-g. 

(Fig. 9) 

To obtain the upper limits on the range of density for this purpose, the largest 3 -hour change in 
the exospheric temperature from the MET model was computed. The exospheric temperatures 
associated with the largest daily 3-hour changes (Ti, T 2 ) are shown in Fig. 10; the range of 
temperature for T E >1000 K is approximately ±200 K. Densities for T E = 1600 K ± 200 K at 400 
and 300 km from Table A are listed in Table 1 . 

(Fig. 10) 
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Table 1. Maximum 3-hour change in thermosphere density for Tr = 1600 ± 200K 


T e (K) 

p 

kg m' 3 ... 

400 km 
Ap 

kg m- 3 

Ap rel 
percent 

p 

kgm' 3 

300 km 
Ap 

kgm' 3 

Ap rel 
percent 

1400 

1.0136 x 10 " 

-0.4118x10'" 

-29 

2.2861x10" 

-2.5213x10'" 

-52 

1600 

1.4324 x 10'" 

0.0 

0 

4.8074 x 10" 

0.0 

0 

1800 

1.8642 x 10 " 

0.4318x10" 

30 

6.9264x10" 

2.1191x10" 

44 


ISS Orbital Decay 

Three different models are used to compute the ISS orbital decay. They are: (Model 1) A 
complete orbital decay model whereby the decay rate is calculated for each orbit by integrating 
the densities (computed using the MET model with the prior 81-day means, prior day value and 
the value of a,, 6. 7-hours prior to the time of interest) at 90 evenly-spaced points around the orbit; 
(Model 2) A simple orbital decay model (described later) using a single average global 
thermosphere density computed from Eq. 7 based on the same inputs to the MET model as 
Model 1 and updated at 3-hour intervals; (Model 3) same as (Model 2) but using the 13-month 
smoothed 10.7 cm solar flux and 13-month smoothed geomagnetic index, ap, in the MET model 
with the daily flux set equal to the 13 -month smoothed solar flux. The results for the complete 
orbital decay model for ISS initial altitudes of 450 and 400 km beginning January 22, 1958, a 
high solar activity condition, are illustrated in Fig. 1 1 . 

(Fig. 11) 

A simple orbital decay model for a circular orbit can be expressed as: 


f * CdB ('. +z )e. v < 


L ^(r a +Z)cos ij 


(14) 
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where r 0 = 6,378,140 m, Z is altitude above the earth, pz is the density at altitude Z (m), V c is the 


satellite circular velocity (m/s) where 


V, = 


l 

B 


L( r o+z)J 


and p = 3.986012 xl0 14 m 3 s' 2 , co = 7.2921 x 10’ 5 rad/s, i = 51.6 deg. The other parameters are 
dz 

as defined in Eq. 11; — is the orbital decay rate (m/s). 


From Model 1 ISS decayed from the initial altitude of 400 km to 253 km in 90 days. From 
Model 2 (Eq. 14), using the same solar inputs as in Model 1, the ISS altitude decayed from 400 
to 263 km in 90 days. From Model 3 the ISS altitude decayed from 400 to 240 km in 90 days 
(lower curves. Fig. 1 1). For the ISS at an initial altitude of 450 km (upper 3 curves in Fig. 1 1) all 
three orbital decay models yield nearly coincident values of decay to 405 km in 90 days. This 
analysis suggests that the simple orbital decay model (Eq. 14 using the average global density) is a 
good approximation for the ISS orbital decay. 

The 100-day ISS orbital decay from initial altitudes of 400 km and 450 km was computed using 
the simple orbital decay model (Eq. 14) with inputs to the MET model of the 81 -day prior 
average solar flux, daily solar flux and ap at 6 hours prior to the time of computation to compute 
the global maximum and minimum exospheric temperature (Eqs. 5 and 6) and then the global 
average density was computed(Eq. 7). The results are shown in Fig. 12 (solid lines) for the years 
1948 to 1994. In a similar manner, the 13 -month smoothed solar flux, and 13 -month smoothed ap 
were used in the MET model to compute the global average density. The daily solar flux was set 
equal to the 13 -month smoothed solar flux. These results for 100 days decay are shown as dots in 
Fig. 12. For initial altitudes of 400 and 450 km there is good agreement between the two MET 
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model inputs(8 1-day prior and 13-month smoothed) to the simple orbital decay model. For solar 
cycle 19 and an initial altitude of 400 km the orbital decay in 100 days is so large that there is a 
risk indicated that the ISS will reenter during high solar activity unless reboosted. In contrast, the 
orbital decay during low solar activity for the initial altitude of 400 km is approximately 10 km in 
100 days. The variation in orbital decay for various initial altitudes and durations is' an important 
consideration in planning re-boost strategies for the ISS. 

(Fig. 12) 

The ISS assembly is to be completed in the year 2002. Assuming that solar cycle 23 began 
near mid-year 1996 places part of the ISS assembly period near the peak of the cycle when solar 
activity effects on orbital decay could be significant. 

CONCLUSIONS 

Based on the MET model for estimation of neutral mass density at ISS orbital altitudes for the 
entire period of record (1947-1994) for the model input variables that produced a calculated 
maximum value for exospheric temperature of 2152.2 K, the altitude above which the 
atmospheric drag acceleration on the ISS would not have exceeded one micro-g is 397 km for 
no uncertainty in the MET model density and 422 km if the MET model underestimates the 
density by 20 percent (the worst case uncertainty) . These altitudes are reduced to 363 km (no 
uncertainty) and 380 km (worst case 20 percent uncertainty) if the maximum exospheric 
temperature is less than 1700 K. The range of this critical altitude versus exospheric temperature 
(Fig. 6) for one micro-g has important implications for the ISS planning strategies. 

For the ISS orbital decay up to 100 days for initial altitudes at 400 and 450 km, the use of the 
13 month smoothed values for the 10.7 cm solar flux and the geomagnetic index (a p ) produces 
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results which are comparable to those produced by the using the specified inputs to the MET 
model. For periods of high solar activity, as in solar cycle 19 (Fig. 12), this study suggests the ISS 
could re-enter the earth’s atmosphere within 100 days from an initial altitude of 400 km if not 
reboosted. 

For development of re-boost operational strategies and for real time operations the complete 
MET model with 81-day running means and daily values for 10.7 cm solar flux and 3-hourly 
values for ap are required. 

The variation in density within an ISS orbit is important information for the assessment of the 
capability of the ISS to meet the micro-g experiment acceleration requirements and requirements 
for torque and attitude control by the control moment gyros. 
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Appendix 

Table A. MET model neutral mass density (kg/m 3 )as a function of altitude and exospheric 
temperature. 


Total Density kg/m 3 
Exospheric Temperature, deg K 

Z (km) 600. 800. 1000. 1200. 1400. 1600. 


250. 

260. 

270. 

280. 

290. 

300. 

310. 

320. 

330. 

340. 

350. 

360. 

370. 

380. 

390. 

400. 

410. 

420. 

430. 

440. 

450. 

460. 

470. 

480. 

490. 

500. 


z ( km) 
250. 
260. 
270. 
280. 
290. 
300. 
310. 
320. 
330. 
340. 
350. 
360. 
370. 
380. 
390. 
400. 
410. 
420. 
430. 
440. 
450. 
460. 
470. 
480. 
490 
500. 


0. 18483E-10 
0 . 1 3158E-1 0 
0. 94599E-11 
0.68587E-11 
0 . 50090E-1 1 
0 . 36812E-1 1 
0 . 27202E-1 1 
0 . 201 99E-1 1 
0 . 1 5064 E-l 1 
0.1127 9E-1 1 
0. 84770E-12 
0 . 6394 5E-12 
0.4841 3E-1 2 
0 . 367 91E-1 2 
0.28069E-12 
0 . 21506E-1 2 
0 . 16555E-12 
0 . 12808E-12 
0. 99660E-13 
0 . 78036E-1 3 
0 . 61 535E-1 3 
0.48899E-13 
0. 39189E-13 
0. 31695E-13 
0. 25886E-13 
0 . 21359E-13 
Total I 


0. 44060E-10 
0. 33150E-10 
0.2517 5E-1 0 
0.19279E-10 
0. 14875E-10 
0. 11553E-10 
0. 90271E-11 
0.70906E-11 
0. 55958E-11 
0.4434 9E-1 1 
0 . 35281E-11 
0 . 28163E-11 
0 . 22551E-11 
0. 18108E-11 
0. 14579E-11 
0. 11766E-11 
0. 95181E-12 
0.77167E-12 
0. 62699E-12 
0 . 51052E-12 
0. 41657E-12 
0 . 34 066E-12 
0 . 27 921E-12 
0 . 22 938E-12 
0 . 18891E-12 
0. 15598E-12 
tensity kg/m 3 


0.7324 8E-1 0 
0.57155E-10 
0.44 972E-1 0 
0 . 35654 E-l 0 
0.284 61E-10 
0 . 22861 E-l 0 
0.184 66E-1 0 
0.14 993E-10 
0. 12230E-10 
0. 10019E-10 
0.82392E-11 
0 . 67 997 E-l 1 
0 . 562 98E-11 
0. 46751E-11 
0 . 38 927 E-l 1 
0. 32494E-11 
0.27186E-11 
0.22794E-11 
0. 19150E-11 
0. 16118E-11 
0. 13589E-11 
0. 11477E-11 
0. 97078E-12 
0. 82241E-12 
0. 69774E-12 
0. 59282E-12 


1800. 

0. 15528E-09 
0 . 13035E-09 
0 . 11028E-09 
0. 93924E-10 
0 . 804 55E-10 
0.69265E-10 
0.59897E-10 
0 . 52000E-10 
0. 45305E-10 
0 . 39598E-10 
0 . 34 7 12E-1 0 
0 . 3051 1 E-l 0 
0 . 26886E-10 
0.23746E-10 
0.21019E-10 
0. 18642E-10 
0. 16567E-10 
0. 14749E-10 
0 . 13153E-10 
0. 1174 9E-10 
0 . 10512E-10 
0. 94181E-11 
0. 84504E-11 
0.75922E-11 
0 . 68300E-11 
0 . 6151 9E-11 


c Temperatu 

2000 . 

0. 16607E-09 
0.14 068E-09 
0 . 12013E-09 
0 . 10330E-09 
0. 89343E-10 
0.77 667 E-10 
0.67816E-10 
0.594 4 5E-10 
0 . 52288E-10 
0. 46135E-10 
0.40820E-10 
0 . 36210E-10 
0 . 321 96E-10 
0.28689E-10 
0.25616E-10 
0.22915E-10 
0 . 20536E-10 
0.184 35E-10 
0. 16575E-10 
0.14 925E-10 
0. 13459E-10 
0 . 12154 E-10 
0. 10989E-10 
0. 99492E-11 
0. 90183E-11 
0 . 8184 0E-11 


e, K 

2200 . 

0. 17454E-09 
0. 14889E-09 
0 . 1281 0E-09 
0 . 11100E-09 
0. 96764E-10 
0.847 94 E-l 0 
0.74637E-10 
0 . 65 953E-1 0 
0.584 7 9E-1 0 
0 . 52009E-10 
0 . 4 6381E-1 0 
0. 41464E-10 
0.37151E-10 
0 . 33356E-10 
0 . 30005E-10 
0 . 27039E-1 0 
0. 24407E-10 
0 . 22065E-10 
0 . 1997 8E-10 
0 . 1 81 1 3E-1 0 
0 . 164 45E-10 
0. 14949E-10 
0. 13606E-10 
0. 12398E-10 
0 . 11311E-10 
0. 10329E-10 


0. 10047E-09 
0. 80505E-10 
0 . 65006E-1 0 
0 . 52856E-10 
0. 43246E-10 
0 . 35584 E-10 
0.294 31E-10 
0.244 57E-10 
0 . 204 12E-10 
0 . 17 1 03E-10 
0 . 14 384E-10 
0 . 121 38E-10 
0 . 1027 5E-10 
0 . 87 22 8E-11 
0 . 7 4 252E-11 
0 . 63365E-11 
0 . 54 200E-11 
0. 46460E-11 
0. 39906E-11 
0 . 34 34 IE-1 1 
0 . 2 9604E-1 1 
0 . 25563E-11 
0 . 221 07E-1 1 
0 . 1 9147E-11 
0 . 16606E-11 
0. 14421E-11 


0 . 12326E-09 
0 . 1007 5E-09 
0. 82965E-10 
0.68775E-10 
0.57348E-10 
0.48074E-10 
0.404 92E-10 
0 . 34 254E-1 0 
0. 29092E-10 
0.247 99E-10 
0 . 21211E-10 
0 . 181 99E-10 
0 . 15660E-10 
0 . 1 3513E-10 
0.1168 9E-10 
0 . 10136E-1 0 
0 . 88081E-11 
0 . 7 67 06E-1 1 
0. 66931E-11 
0. 58509E-11 
0.51237E-11 
0.44 941E-1 1 
0 . 394 80E-11 
0 . 34733E-1 1 
0 . 30600E-11 
0 . 26993E-1 1 


0.14131E-09 
0.11726E-09 
0. 98035E-10 
0 . 82503E-1 0 
0 . 6 9832E-1 0 
0.59409E-10 
0 . 5077 IE-1 0 
0.43567E-10 
0.37524E-10 
0. 32429E-10 
0.2811 4 E-l 0 
0.2444 4 E-l 0 
0. 21309E-10 
0. 18624E-10 
0 . 1 6315E-1 0 
0.14 324 E-l 0 
0 . 12602E-10 
0 . 11109E-10 
0. 98110E-11 
0.867 98E-11 
0 . 7 6918E-1 1 
0 . 68270E-1 1 
0 . 60685E-11 
0.54019E-11 
0 . 4 8150E-1 1 
0.4297 3E-1 1 
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